USTC-ICTS-08-23 



Reconstructing f(R) Theory from Ricci Dark 

Energy 

Chao-Jun Feng 

Institute of Theoretical Physics, CAS, 
Beijing 100080, P.R.China 
Interdisciplinary Center for Theoretical Study, USTC, 
Hefei, Anhui 230026, P.R.China 

fengcj@itp.ac.cn 



In this letter, we regard the f(R) theory as an effective description for the 
acceleration of the universe and reconstruct the function f(R) from the Ricci 
dark energy, which respects holographic principle of quantum gravity. By 
using different parameter a in RDE, we show the behaviors of reconstructed 
f(R) and find they are much different in the future. 



1 Introduction 



The accelerating cosmic expansion first inferred from the observations of 
distant type la supernovae [1] has been strongly confirmed by some other in- 
dependent observations, such as the cosmic microwave background radiation 
(CMBR) [2] and Sloan Digital Sky Survey (SDSS) [3]. An exotic form of neg- 
ative pressure matter called dark energy is used to explain this acceleration. 
The simplest candidate of dark energy is the cosmological constant A, whose 
energy density remains constant with time p\ = A/8irG and whose equation 
of motion is also fixed, = Pa/ Pa = — 1 (-Pa is the pressure) during the 
evolution of the universe. The cosmological model that consists of a mixture 
of the cosmological constant and cold dark matter is called LCDM model, 
which provides an excellent explanation for the acceleration of the universe 
phenomenon and other existing observational data. However, as is well know, 
this model faces two difficulties, namely, the 'fine-tuning' problem and the 
'cosmic coincidence' problem. The former also states: Why the cosmological 
constant observed today is so much smaller than the Plank scale, while the 
latter states: Since the energy densities of dark energy and dark matter scale 
so differently during the expansion of the universe, why they are at the same 
order today? 

To alleviate or even solve these two problems, many dynamic dark en- 
ergy models were proposed such as the quintessence model relying on a scalar 
field minimally interacting with Einstein gravity. Here 'dynamic' means that 
the equation of state of the dark energy is no longer a constant but slightly 
evolves with time. So far a wide variety of scalar-field dark energy models has 
been proposed including quintessence mentioned above, k-essence, tachyons, 
phantoms, ghost condensates and quintom etc. .Despite considerable works 
on understanding the dark energy have been done, we can not answer this 
question at present, because we do not entirely understand the nature of dark 
energy before a complete theory of quantum gravity is established, since the 
dark energy problem may be in principle a problem belonging to quantum 
gravity [1]. Actually, there is still a different way to face the problem of cos- 
mic acceleration. Since general relativity is only tested within solar system 
up to now, it is possible that the observed acceleration is not the manifesta- 
tion of another ingredient in the cosmic pie, but rather the first signal of a 
breakdown of our understanding of the laws of gravitation, as stressed by Lue 
et al.[5|. From this point of view, one may consider the modification to the 
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Einstein-Hilbert action at larger scales with higher order curvature invariant 
terms such as R 2 , R^ U R^ U , R^^RpvaPi or R\3 k R as well as nonminimally 
coupled terms between scalar fields and geometry (such as <j) 2 R). These terms 
naturally emerge as quantum corrections in the low energy effective action 
of quantum gravity or string theory [3 [6]. The interesting models following 
this line include f(R) and DGP gravity, and in this letter, we will focus on 
f(R) theory where the modification is a function of the Ricci scalar only 

Although we are lacking a quantum gravity theory today, we can still 
make some attempts to probe the nature of dark energy according to some 
principle of quantum gravity [7J. It is well known that the holographic prin- 
ciple is an important result of the recent researches for exploring the quan- 
tum gravity(or string theory) j3] . So that the holographic dark energy model 
(HDE) constructed in light of the holographic principle possesses some sig- 
nificant features of an underlying theory of dark energy [7J. Recently, Gao 
et.al [9] proposed a holographic dark energy model in which the future event 
horizon is replaced by the inverse of the Ricci scalar curvature, and they call 
this model the Ricci dark energy model(RDE). Of course, this model also 
respect the holographic principle. 

In this letter, we regard RDE as the underlying theory of dark energy 
and reconstruct the corresponding f(R) theory as an equivalent description 
without resorting to any additional dark energy component, namely RDE 
is effectively described by f(R) theory. In Section II, we will briefly review 
RDE and f(R) models, and reconstruct function f(R) from RDE model in 
Section III. In the last section we will give some conclusions. 



2 Briefly Review on RDE and f(R) theory 

Holographic principle [10] regards black holes as the maximally entropic ob- 
jects of a given region and postulates that the maximum entropy inside this 
region behaves non-extensively, growing only as its surface area. Hence the 
number of independent degrees of freedom is bounded by the surface area in 
Planck units, so an effective field theory with UV cutoff A in a box of size L 
is not self consistent, if it does not satisfy the Bekenstein entropy bound [11] 
(LA) 3 < Sbh = ^L 2 M^, where MZ 2 = G is the Planck mass and Sbh is the 
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entropy of a black hole of radius L which acts as an IR cutoff. Cohen et.al. 
[T2"] suggested that the total energy in a region of size L should not exceed 
the mass of a black hole of the same size, namely L 3 A 4 < LM 2 . Therefore 

the maximum entropy is S^h- Under this assumption, Li [13] proposed the 
holographic dark energy as follows 

p A = 3c 2 M p 2 L- 2 (1) 

where c 2 is a dimensionless constant. Since the holographic dark energy with 
Hubble horizon as its IR cutoff does not give an accelerating universe [H] , Li 
suggested to use the future event horizon instead of Hubble horizon and par- 
ticle horizon, then this model gives an accelerating universe and is consistent 
with current observation [T3"t [T5] . For the recent works on holographic dark 
energy, see ref. [IE]. In the following, we are using units 8nG = c = K = 1. 

Recently, Gao et.al [9] proposed a holographic dark energy model in which 
the future event horizon is replaced by the inverse of the Ricci scalar curva- 
ture, and they call this model the Ricci dark energy model(RDE). This model 
does not only avoid the causality problem and is phenomenologically viable, 
but also solve the coincidence problem of dark energy. The Ricci curvature 
of FRW universe is given by 

R=-6(H + 2H 2 + ^-), (2) 



a 2 

where dot denotes a derivative with respect to time t and k is the spatial 
curvature. They introduced a holographic dark energy proportional to the 
Ricci scalar 

Px = 3aUl + 2H 2 + ^\ oc R (3) 

where the dimensionless coefficient a will be determined by observations and 
they call this model the Ricci dark energy model. Solving the Friedmann 
equation they find the result 



PX <* n -3x , (4-2)3 (A s 

where Q m o = p m o/3-ffo, x = In a and fo is an integration constant. Substi- 
tuting the expression of px into the conservation equation of energy, 

1 dp x ( , 

PX = ~ PX ~ 3^7 (5) 
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we get the pressure of dark energy 



Px = 



{fa ~ I) ^ 



(6) 



Taking the observation values of parameters they find the a ~ 0.46 and 
/o ~ 0.65 [9j. The evolution of the equation of state wx = Px/px of 
dark energy is the following. At high redshifts the value of wx is closed 
to zero, namely the dark energy behaves like the cold dark matter, and 
nowadays wx approaches —1 as required and in the future the dark energy 
will be phantom. The energy density of RDE during big bang nucleosynthe- 
sis(BBN) is so much smaller than that of other components of the universe 
[Qx\iMeV < 10~ 6 <C 0.1 when a < 1), so it does not affect BBN procedure. 
Further more this model can avoid the age problem and the causality prob- 
lem. 

Consider the modification of gravity, one can add terms like R 2 , R^R^, 
RV VOi PR pva p } R\3 k R or nonminimally coupled terms to the effective Lagrangian 
of gravitational field when quantum corrections are considered. In f(R) the- 
ory, the modification is adding a function of the Ricci scalar only, and the 
action of it is as follows 



where C m is the matter Lagrangian, and f(R) is a function of R. Then 
we want to obtain the modified Friedmann equations by varying the gen- 
eralized Lagrangian. However, it is not clear how the variation has to be 
performed [18]. Assuming the FRW metric, the equations governing the dy- 
namics of the universe are different depending on whether one varies with 
respect to the metric only or with respect to the metric and the connection. 
These two possibilities are usually called the metric and the Palatini [19J ap- 
proach respectively. It is only in the case of Einstein gravity f(R) = R 
that these two methods give the same result. The problem of which method 
should be used is still a open question and a definitive answer is likely far to 
come. In ref. [18], a method was proposed to reconstruct the form of f(R) 
from a given Hubble parameter H(z) from observational data such as SN's 
Gold datapQ in the metric formulation. What is needed to reconstruct f(R) 
in their approach is an expression for H(z), so we can use a H(z) predicted by 
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a given dark energy model to determine what is the f(R) theory which give 
rise to the same dynamics[l8j- In the following, we will follow this method 
to reconstruct f(R) from RDE, so we use the metric formulation, although 
the dynamical equations are more simply in Palatini formulation. 



Variation with respect to the metric leads to the modified Einstein equa- 
tion [2D] 

G^u = Rpuv — —g^R = T^ u ^ + Tp V , (8) 
where is the Einstein tensor and 

T^ rv) = lf(R) - Rf'(R)} + f'(RY aP (9»«9u P - g^gaf,)} (9) 

and the stress-energy tensor of matter 

= f^/fXR) (10) 

with T™ the standard minimally coupled matter stress-energy tensor. Here 
and in the following, we denote with a prime the derivative with respect to 
R. With the FRW metric, we obtain the modified Friedmann equations 



tt2 k 1 

a 2 3 



Pin 

Pcurv + 



(11) 



k 

2^- + H 2 + — = - ( Pcurv + Pm ) , (12) 



a a 



where p m and p m are the matter-energy density and pressure respectively, 
and we have defined the same quantities for the effective curvature fluid as: 

Pcurv = — r^T { \ [/(^ - Rf'(R)} ~ SHRf'(R) | (13) 



and 



Pcur V = j^j [\Rf"(R) + Rf"(R) + R 2 f"\R) - \ lf(R) - Rf'(R)]] ■ 

(14) 

Applying the Bianchi identity to eq.flU), we obtain the conservation law for 
the total energy density p to t = Pcurv + Pm/ f'(R) as follows 

ptot + 3H(p t ot + Ptot) = . (15) 
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In the following, we will study the case of flat universe, i.e. the spatial 
curvature k = with the matter as dust, namely p m = and we do not 
consider the interaction between the matter and the curvature fluid. Thus, 
the matter energy density is conserved so that p m = 3HQfl mQ e~ 3x . 

In fact, eq. ffTTj) . f|T2l and ffTBT) are not independent, so we will consider 
only eq. ffTTj) and ffTo^) . Combine these two equations we obtain 

H = -^(^{zHfcrnoe- 3 * + Rf'(R) + R [Rf'"(R) - Hf'(R)} j . (16) 

Using the relation d/dt = Hd/dx to replace the variable t by x, eq. (1161) can 
be rewritten as a third order differential equation of f(x) 

C 3 (a;)0 + C 2 (x)0 + C x (x)£ = -3O m0 e- 3 * , (17) 

where C n (x) consists of h(x) = H(x)/H and its derivatives, where H is 
the present Hubble parameter,see Appendix A. From eq. fll7l) . one can see 
that what is needed to reconstruct f{R) is an expression for h(x). As a 
consequence, one could adopt for h(x) predicted by a given dark energy 
model and determine what is the f{R) theory [IS]. In the following, we will 
reconstruct f{R) according to RDE. 

3 Reconstructing f(R) from RDE 

From eq.(8) in [9], we get 

h\x) = -^n m0 e- 3x + foe-^) x , (18) 
1 — a 

in RDE with h(x = 0) = 1 as definition. Thus 

1 On0 + /0 = l. (19) 



2 - a 



The differential equation (ITTI) is so complicated that we will solve it nu- 
merically. According to ref.JTS], the boundary conditions, i.e. the values of 



7 



/ and its first and second derivatives with respect to x evaluated at x = 
are as follows 



dx , ;: .,, 



11 

f(x = 0) = f(R ) 



dR 

dx 

(PR 

dx 2 



2 = 



2=0 



6H 2 



m0, 



Rn 



(20) 

(21) 
(22) 



These conditions are chosen on the basis of physical consideration only. 
Rewrite eg. (11 II) explicitly with 87rG and k = as 



H 2 



8ttG 



Pc 



+ 



Pr. 



f'(R) 



(23) 



This equation shows that the function f'(R) is equivalent to redefine the 
Newton gravitational constant G as Gj f'(R), that is time dependent in f(R) 
theory. In order to be consistent with solar system experiments at x — 0, the 
effective gravitational constant G/ f'(Ro) must equal to G, thus f'(Ro) = 1, 
so we get 



f'(Ro 



dR 

dx 



dx 



(24) 



2 = 



which leads to eq. f)20p . A second condition comes from that any f(R) theory 
must fulfill the condition f"(Ro) = in order to not contradict solar system 
tests [2TJ, then it gives rise to eq. (l2Tl) . Finally, the present value of p curv in 
eq.flHP is 

p curv (x = 0) = f{R °: I' R ° . (25) 



By using the eq. ffTT]) evaluated at present (x 
the final boundary condition eq. (!22|) . 



0), and eq. tpB]) . we obtain 



Given a = 0.46, Q m o = 0.27 and h(x) in eq. ([T8|) . with relation eq. ([T9!) 
and three boundary condition eq. (l2~0l) - (l2~2l) . the differential equation eq.( fl7l) 
can be solved numerically. We plot the function f(R) with respect to R in 
Fig.l. 



8 



Reconstruct f(R) 
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Figure 1. Reconstructed f(R) with < z < 10, where the redshift 
z = e~ x — 1, Q m0 = 0.27 and a = 0.46 (solid), a = 0.5 (dashed) and 
a = 0.54 (dash- dotted). 

From Fig.l one can see that for small \R\ (small z also), the functions are 
distinguishable for different parameter a. Differences between these function 
f{R) become significant when \R\ (or z) increases. In order to compare with 
the results in [IB] , we also show our results on a If — IR plane in Fig.2, where 
If = ln(— /) and IR = ln(— R) used in [TB], and our results are consistent 
with theirs. Moreover, Fig.l and Fig.2 indicate the parameter a plays a im- 
portant role in the remote past. Actually, the value of a also determines the 
future evolution of R. To illustrate this, we plot the evolution of R in the 
future in Fig.3. As is expected, for a < 0.5, the curves indicate \R\ — > oo in 
the future, which is the behavior of phantom with equation of state smaller 
than —1 dominating over others. The energy density of phantom increases 
with time, tears apart structures and a Big Rip is unavoidable. For a = 0.5, 
\R\ varies a little, and the dark energy becomes more and more like a cos- 
mological constant. For a > 0.5, R will vanish in the future. 
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Reconstruct lf-lR 
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Figure 2. Reconstructed f(R) in If — IR plane with < z < 10, where the 
redshift z = e~ x — 1, Q m o = 0.27 and a = 0.46 (solid), a = 0.5 (dashed) and 
a = 0.54 (dash-dotted). If = ln(-/) and IR = ln(-R). 

R(z) 
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Figure 3. The future evolution of R with < z < 10, where the redshift 
z = e~ x — 1, Q m o = 0.27 and a = 0.46 (solid), a = 0.5 (dashed) and 
a = 0.54 (dash-dotted). 
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Reconstruct f(R) a=0.46 




Figure 4. Reconstructed f(R) with the redshift z = e~ x — 1 from around 2 
down to —1, Q m0 = 0.27 and a = 0.46. The arrow denotes the decreasing 
direction of z, and the point corresponds to the current value at z = 0. 



Reconstruct f(R) a=0.5 
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Figure 5. Reconstructed f(R) with the redshift z = e~ x — 1 from around 2 
down to —1, Q m0 = 0.27 and a = 0.5. The arrow denotes the decreasing 
direction of z, and the point corresponds to the current value at z = 0. 
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Figure 6. Reconstructed f(R) with the redshift z = e~ x — 1 from around 2 
down to —1, Q m o = 0.27 and a = 0.54. The arrow denotes the decreasing 
direction of z, and the point corresponds to the current value at z = 0. 

In Fig. 4, Fig. 5 and Fig. 6, we can see that the difference is more distinc- 
tively reflected by the function f(R) reconstructed according to the future 
evolution of RDE. For a = 0.46, there exists a turnaround point on the 
curve at present epoch (z = 0), and the decreasing \R\ begin to increase 
at this point, which means the phantom-like dark energy will dominate the 
universe. The existence of the turnaround point is a common feature for 
all the phantom-dark energy models realized in f(R) theory because of the 
competition between dark energy and matter [22] . For a = 0.5, f(R) linearly 
depends on R up to a constant, which corresponding to the de Sitter space, 
where f(R) = R + 2A. For a = 0.54, f(R) increases from negative to positive 
and seems a inverse power law dependence on R in the future. 

4 Conclusions 

In conclusion, we have followed the method proposed in ref.|18j to recon- 
struct the function f(R) in the extended theory of gravity according to 
the Ricci dark energy model, which respects holographic principle of quan- 
tum gravity. We show the behaviors of f(R) reconstructed with parameter 
a = 0.46, 0.5, 0.54 in RDE and find that the dependence of f(R) on R is 
different for different a, and such a difference is much more distinctive in the 



Reconstruct f(R) o'=0.54 
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future. The basic reconstruction procedure is simply: once the function h(x) 
given by some dark energy model, one can solve the differential equation (1171) 
to obtain / with relation eq. ffT9l) and three boundary condition eq. ([20]) -( T22|) . 
In our case, h(x) is given by RDE in ffl8|) and results of reconstruction are 
shown in Fig. 1-6. The parameter a plays a important role to determine the 
dependence of / on R, so we hope that the future high precision observation 
data may be able to determine it and reveal some significant features of the 
underlying theory of dark energy. 

It should be noted that RDE is obtained within the framework of general 
relativity, rather than any other extended gravity theory such as f(R) theory. 
What we done in this letter is to reconstruct the f(R) theory to effectively 
describe RDE in Einstein gravity. Whether RDE can be generalized to f(R) 
theories is question worth further investigation as HDE [22J. 
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Appendix A 



The coefficients C n {x) in eq. ([T71) are 



d = 2h 2 



d 2 R 

dx 2 



dR\ 
dx J 



h 



, d 3 R 

dx 3 



+ 



ldh 2 
2 dx 



ti; 



d 2 R 

dx 2 



dR s 

dx 



C 2 



-2h 2 



d 2 R 

dx 2 



^ dh 2 (dR^ 



dx V dx 



fdR\ 
\dx J 



-2 



ldh 2 



2 dx 



dR\ 
dx J 



-i 



c. 



\dx J 



(26) 
(27) 

(28) 
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R(x) in eq.© in the case of k = is 

2 j- 

2 rfa; 
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